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1 Introduction and description of results 

The Airy kernel 

Ai(x)Ai'(y) - Ai(;/)Ai'(a:) 

KAi(a;,y) :== 

x-y 

has been known since the nineties (see [28], [13] and [29]). It was introduced in relation with 

the study of random matrices. Its determinant, in fact, has been used to describe the statistical 

behavior of the Gaussian Unitary Ensemble near the edge of the spectrum. One of its most striking 

features has been discovered already in 1991 by Tracy and Widom [31]. Consider, for simplicity, the 



operator related to the Airy kernel restricted to the semi-interval [s, oo). Its Frcdholm determinant, 
depending on the point s, can be expressed through the formula 

det fid- (i\:Ai)|[s,oo)j =exp - / {x - s)q{x)^dx\ (1.1) 

where q{x) is the Hasting-McLeod solution of the Painlevc II equation 

q" = 2q^ + xq 

uniquely determined by its asymptotics at infinity q(x) ^ A[{x) (for an alternative proof of this 
result, based on the theory of solitonic equations, see also [3]). On the other hand , after its in- 
troduction in the nineties, mathematicians and physicists discovered that the Airy kernel and the 
Tracy- Widom distribution (i.e. its determinant) are not related just to GUE but, rather, they are 
a sort of universal objects appearing in many different statistical models such as, just to cite few 
of them, one-dimensional non intersecting Brownian motions, random growth models and random 
partitions (see for instance [32], [26], [12], [5], [17]) 

Few years later it appeared that also the Pearcey kernel has similar properties. It was introduced, 
around 1998, in relation with matrix models with external source [14], [10], 1-dimensional Brownian 
motions [10] and plane partitions [30]. Tracy and Widom, in [33], deduced some differential equa- 
tions for its determinant and, with a different approach, some other PDEs have been obtained in 
[4] and [2]. In particular, in [2], a general method to find some PDEs for general (l,p)-kerncls is 
described. The Airy and the Pearcey kernel correspond to the case p ~ 2,3 while for p > 3 it is 
believed, as reported in [7], that these (l,p)-kernels describe, after a proper rescaling, the statistical 
behavior of multi-matrix models near the edge of the spectrum. These kernels are of the form 

Iff ge,(M)-e„(A) 

'""^'■'•>-i2wL"''L'"^^^7r- C'^' 

acting on / := [oi, 02] U . . . U [a^-i, o-n]] here the contours of integration 71, 72 are curves extending 
to infinity or formal linear combinations thereof, and they are chosen so that the real part of 
the exponential tends to minus infinity and hence yields absolutely convergent integrals. It is to 
be remarked -however- that no independent proof of the total positivity of these kernels exists 
in the literature and hence it is not clear whether they describe a determinantal point process. 
Nevertheless the treatment of their Fredholm determinants can be addressed by completely similar 
methods and hence we find it worthwhile mentioning them here. 



Fredholm determinants of the type (1.1) with Ka\ replaced by K'^p^ (1-2) are instances of de- 
terminants of operators with integrable kernels in the sense of Its-Izergin-Korepin-Slavnov [22], as 
we recall presently (a very concise and clear account is given in [20]). 

Given a piecewise smooth oriented curve C on the complex plane (possibly extending to infinity) 
and two p- vectors /(A),^(A) £ LP'iC) we define the kernel K as 

A — /i 

Since the kernels we will treat will always satisfy this condition we assume that C has no self 
intersections. We say that such kernel is integrable if f^{\)g{\) = (so that it is non-singular on 
the diagonal). We are interested in the operator K : L?{C) — J- L^{C) acting on functions via the 
formula 



{Kf){\)= j KiX,fi)f{fi)dfi 

and, in particular, we are interested in the Fredholm determinant det(Id — iiT). The key observation 
is that, denoting with d the differentiation with respect to any auxiliary parameter on which K 
may depend, we obtain the formula 

aiogdet(Id-ii') = -Triild + R)dK) (1.4) 

where R is the resolvent operator, defined as i? = (Id — K)^^K. Moreover R is again an integrable 
operator, i.e. 

A /i 

and F, G can be found solving the following RH problem: 

r+(A)=r_(A)A/(A) Xec (1.5) 

T{X)=1 + 0{X-^) A ^00 (1.6) 

M(A)= l-27nf(X)f{X) (1.7) 

F(A)= r(A)/(A) G(A) = (r-i)T(A).9(A). (1.8) 

In several cases of interest the Fredholm determinant for such a kernel coincides with the notion 
of isomonodromic tau function introduced in the eighties by Jimbo Miwa and Ueno [25, 23, 24] to 
study monodromy-preserving deformations of rational connections on P^. In particular in [11] the 
authors proved that this is the case'* for a large class of kernels including the so-called 2-F'i-kernel, 



■^More precisely Borodin and Deift proved that the Fredholm determinants they considered are isomonodromic 
tau functions for a certain system of Schlesinger equations. 



the Jacobi and the Whittaker kernel Another good host of examples is provided in [20] and -to a 
certain extent- Airy and Pearcey and the other to be considered here have a large overlap with the 
situation addressed ibidem. 

In the paper we will first connect the Fredholm determinants of these integrable kernels in general 
terms with the tau function associated to Riemann Hilbert problems introduced in [8] following [27] , 
which reduces to the Jimbo-Miwa-Ueno definition of isomonodromic tau function when the relevant 
RHP comes from the (generalized) monodromy problem of a rational ODE; this is achieved in 
Sect. 2 and in particular in Thm. 2.1. The approach followed in the previous literature to connect 
Fredholm determinants of the Airy kernel on a collection of intervals / = [ai, 02] U . . . U [ajv-i, ojv] 
was to set-up a RH problem with jumps on / for instance in [19], to study the integrable differential 
equations related to the Airy kernel from the point of view of isomonodromic deformations. We 
take a different approach and define a new {N + 1) x {N + 1) RH problem with jumps on the 
contour C :~ 71 U 72. This RH problem is related to a certain operator acting on L^{C) and, 
through an appropriate Fourier transform, we prove that such operator has the same determinant 
as the operators related to the (l,p)-kernel. This RHP is quite suitable for a Deift-Zhou steepest 
descent analysis in asymptotic regimes. As applications of our approach we give one more (quite 
straightforward) proof of the Tracy- Widom result for the Airy kernel, we study the asymptotics of 
the Pearcey kernel using the non-linear steepest descent method and we find a Lax Pair for the 
Pearcey process. 

The Lax formalism is developed in Sect. 4.2; using the isomonodromic method [25] we re-derive 
the recent nonlinear PDEs satisfied by the Fredholm determinant of the Pearcey kernel for one 
interval [2]. We find also a new and independent PDE (as to be expected) in Prop. 4.2. While the 
equation itself is rather complicated, the conceptual significance is that the general solution of the 
three equations depends on a finite number of parameters, rather than functional ones. 

The second main goal of this paper is to relate in a precise form the asymptotic of the Pearcey 
process and see how it "becomes" an Airy process; more precisely, using our setup we consider the 
asymptotic behavior of the Fredholm determinant of the Pearcey kernel when the endpoints of the 
intervals is very large as the time parameter grows; we show that the determinant is asymptotically 
factorized into the product of Fredholm determinants of the Airy process. While this is to be 
expected on the ground of physical considerations, we believe this is the first mathematically rigorous 
proof of this factorization; a similar but simpler case has been previously studied in [1] using the 
ordinary steepest descent method instead of the non-linear one used here. 

The proof is completely detailed for the case of a single large interval in Theorem 5.1, while 
the case of several intervals is stated in Theorem 5.2 and is not significantly different (as it will be 
apparent) and hence left to the reader. 



2 Fredholm determinants of integrable kernels and tau func- 
tions of Riemann— Hilbert problems 

Following [8] one may consider a slightly more general notion of tau-function associated to any 
Riemann-Hilbert problem (RHP) depending on parameters and which reduces to that of Jimbo- 
Miwa-Ueno in case such a RHP coincides with the one associated to a rational ODE. We briefly 
recall the setup of [8] ; suppose a RHP is posed on a collection of oriented contours C 

r+{X;s}^r^iX;s)M{X;s) (2.1) 

r{X;s}^l + 0{X-^) , A^oo. (2.2) 

where M{X; s) : C ^>- SLr(C) are some suitably smooth functions of A, depending smoothly (analyt- 
ically) on additional deformation parameters, denoted here generically by s. On the space of these 
deformation parameters, we introduce the following one-form^ (here and below we will denote with 
' the derivative with respect to A) 

c.M(a):=^Tr(rii(A)r(A)Sa(A))^ (2.3) 

Eg{X):^dM{X)M-\X) (2.4) 

The definition (2.3) is posed for arbitrary jump matrices; in the case of the RHP (1.5)-(1.6)-(1.7) 
the spontaneous question arises as to whether lum in (2.3) and the Fredholm determinant are related. 
The answer is positive within a certain explicit correction term, see Thm. 2.1. 

In [8] it was shown that ljm is also the logarithmic total differential of the isomonodromic tau 
function of Jimbo-Miwa-Ueno (in the cases of RHPs that correspond to rational ODEs). 



Theorem 2.1 Let 


f{X;s),g{X;s) : C; 


K S ^ C 


and consider the RHP with 


jumps as 


m (1.7). 


Given 


any vector field d in the space 


of the parameters S of thi 


z integrable 


kernel we 


have the 


equality 




















wm(9) == 


a In dot (Id 


-K)- 


-H(M) 






(2.5) 


where 


'jJAiid) is as 


in (2.3) and 
















H{M) : 


= Hi{M)-H2{M) 


-IM^ 


'g + f'^dgyX- 


2m 1 ff'dffdX 




where 


' denotes differentiation with respect to X. 













^In [8] the sign is the opposite, which would lead to the tau function to have poles instead of zeroes. We correct 
this here. 



The proof is found in Appendix A. In the case we will treat in the following sections, moreover, 
we will have H(M) ~ 0. Hence it is possible to define, up to normalization, the isomonodromic tau 
function tjmu '■— exp(/aJM) and this object, thanks to the previous theorem, will coincide with 
the Fredholm determinant det(Id — K). 



3 The Airy kernel 

We start considering the Airy kernel 
KAi{x,y) 



dfi / dA- 



^i)Ap)-Sy{\) 



X 



^x(A^) 



3 



Xfl 



(3.1) 



(27ri)2 j^^ j^^ A-fj, 

where jr is a contour in the right half-plane which extends to infinity along the rays arg(A) = 
±^ and 7l = ^7_r (the contours are as in Figure 1). We consider the Fredholm determinant 
det(Id — Kj^iXi) where we denote with the same symbol the Airy kernel and the related operator 
acting on L^(M), xi is the characteristic function of the collection of intervals / :— [ai, a2]U [03, a4]U 
. . . U [flAT-i, a^] for even N and / 
remark is self evident: 



[ai, 02] U [03, 04] U . . . U [ajv, 00) for odd N. The following 



Remark 3.1 Let's denote with Ka{x,y) — KAi{x,y)xia,oo){x). Then we have that 



N 

KUx,y)xi{x) = Y.^-iy^'K,^{x,y). 






(3.2) 



Our goal is to setup a RHP associated to the Fredholm determinant of K{x,y)xi{y)', of course 
the interest is somewhat limited since this issue has been thoroughly investigated in the literature 
(see for instance [15] for the similar case of the sine kernel and more generally [11]). We point out 
-however- that the RHP that we are going to set up is not of the same nature as the natural one 
considered in [19] or in [11]. We will use this formulation later on to investigate the asymptotic 
behavior of the Pearcey Fredholm determinant. 

Definition 3.1 Given I = [01,02] U . . . U [ajv-i,ajv] for even N or I — [01,02] U . . . U [aAr,cx)) for 
odd N we define the related AiO-RH problem (standing for Airy-operator Riemann-Hilbert problem) 

r+(A)-r_(A)(l-G(A)) Ae7:=7;^U7L 



r(A) = l+C'(A-i) A 



G(A) 



e''"i(^) -e^.^2W 




(-) 



A^+lgi?o 



rW 



X7rW 



C-^.i(A) 

c-^a.,W 



-^a^W 



(3.3) 



X-TL W 



which consists of finding an analytic matrix-valued function T{X) on C\"f normalized to the identity 
at infinity and such that the limiting value ofT{X) approaching the contour from the left and from 
the right (T_|_(A) and r_(A) respectively) are related through the jump matrix G{X). 

Remark 3.2 (Solvability of the problem in Def. 3.1) The problem in Def 3.1 will be asso- 
ciated to an integrable kernel in Prop. 3.1. The solvability is then equivalent (by the theory of 
Its-Izergin-Korepin-Slavnov [22]) to the non-vanishing of the corresponding Fredholm determinant. 
It will be shown in Thm. 3.1 that this determinant is precisely the gap-probability of the Airy 
process, hence (strictly) positive. This guarantees automatically the solvability for real values of 
the parameters a.j. In the case N ^ 1 the solvability of the problem in Def. 3.1 is equivalent to 
the absence of poles in the Hastings-McLeod solution, shown in [21]. The sam,e remark, mutatis 
mutandis, shall apply to the Riemann-Hilbert problem in Def. 4-1 using Theorem ^.1. 

The RHP of Def 3.1 is associated to a certain Fredliolm determinant of an integrable operator 
on 7fl U 7l as explained in the introduction. 

Proposition 3.1 The Riemann-Hilbert problem (3.3) is the RH problem associated to the integrable 
operator 



K{X,fi) 



/(A) 



A — /I 





/ 


r ^^ 1 

eir 


1 







2TTi 








\ 






X7«(A) + 





p — 6-4-aiA 



e~ -Q--^o.N^ 



\ 



X-tlW 



(3.4) 



(3.5) 



.9(A) 



g-g-— aiA 



(-1) 



N+l^^-a^X 



X7h(A) 



e s 




X7l(A) 



(3.6) 



Proof: It is just enough to verify that l-27r«/(A)5(A)'^ = l-G(A). Q.E.D 

Remark 3.3 The precise form of the contours jl, jr is not essential; in fact we could replace 7^ 
by the imaginary axis. One can directly show that the Riemann-Hilbert problems obtained by this 
"contour deformation" of 7^ to iR are equivalent. 

Theorem 3.1 The following identity holds 

det(Id - K) = det(Id - Kmxi)- 

where the operator defined by K is a trace-class operator on L'^{'^l U 77?). 



7 



Proof: We start with observing that any operator on L'^{'^l U 7_r) — LF'{il) ffi L'^{ir) = "Hi ® 7^2 
can be written as a 2 x 2 matrix of operators with {i,j) entry given by an operator "Hi — > "Hj. 
Writing out K in fuU we have 



K{\y) 



2i7r A — /i 



Define operators J^, (?„ as follows 

5(A) 



H^ 



H^ 



si 1 

6 

2m 



LHil) 



e^-./Md^ 



7h 



C-M 



c^ / ^-Ai g(A) 
2m 



e 6 



71. 



fl — \ 



dX. 



(3.7) 



(3.8) 



(3.9) 



We will think of these operator as acting on H := Hi (B H2 by extending them trivially to the 
orthogonal complements of the respective domains. The operators Qa and T are of Hilbert-Schmidt 
class (HS for short) in "H: this follows from the convergence of the following expressions: 



l^a||-«,2 






< +00 



(3.10) 



m 



H.2 



j'''ir^(-.i 



-) 



< +00 



(3.11) 



Remark 3.4 Note that they would still be of HS class even if we replaced 7l by iR (which we will 
do later). 

More is true: both Qa and T arc of trace-class. To see this it is sufficient to write them as 
the composition of two HS operators. To achieve this goal we introduce an additional contour 
7o := zM + e not intersecting cither of 7L^fl: and extend the Hilbert space by adding H :— L'^{jo)(BH. 
Consider now the two operators on H defined by 

(3.12) 



ci'\fm 



/(^)cTi^-"^d;U 



2^^(C - m) 



C?n/i)(0:=c-ir+< 



fe(C)dC 
2tTT{C - 



and extended trivially on the orthogonal complements of the respective domains within H. They 
both are HS in Ti because 



cl^) 



W,2 



70 



g25R(V-aM) 

IdCI / IdMl 4^^^^^^ < +00 



7-R. 



cp) 



W,2 



7-L 






70 



A simple Cauchy residue-computation closing the integration by a big circle to the left shows 

GafiO=Ci'Kci'^f{0 (3.13) 



and hence Qa is the composition of two HS operators and thus of trace-class. For T one has to use 
instead the operators 

r hjQdc 

/70 2*^(C - m) 



VW : L^jl) ^ i"(7o) 
2?(^H/)(C) := 



e-^ /(A)dA 



2z^(A-C) 
(both HS) that realize J" = X>(2) „ ^(i) 



p(2)(/^)(^) :^e'^ 



(3.14) 



Remark 3.5 Note that if we replace jl by iR it is stiU true that Ga = C^/^ o ci^\ T = P^^^ o X'(i) 
but now 2?'^^ and Co /aiZ to 6e ifiS, and hence we do not know whether Ga,^ <ire of trace-class 
between L'^{jjf) and L^(ilR) (they are still HS, see Rem. 3.4)- 

According to the split H = L'^{'^l) ® L/^ilR) and using matrix notation, we can write det(Id — K) 



det 



Id- 



ELi(-P+'5a, 



J^ 







N 



det[Id-^(-l)-'+igo^. oj-] 



(3.15) 



Both determinants are well-defined because of the form Id+trace class. This identity between 
Fredholm determinants of operators follows by multiplying on the left the operator in the left hand 
side by the operator (with unit determinant) 

(3.16) 



Id 











Z^7 = l( ) £^aj 







Note that the operator /C := X]i=i(^l)"'^^^ai ° -^ appearing in the second term of (3.15) is an 
operator acting on ^^(7^) into itself and with kernel 



N 



/C(e,A):=5](-l)^+i- 



Vaji 



3 = 1 



(2i7r) 



d/i 



Q 3 u,j^ 



IR 



(e-M)(M-A) 



A,ee7L 



(3.17) 



The integral operator on L'^{jl) defined by JC can be made to act on L^(iR), with the same kernel 
but where now ^, A G iM: let us denote this new operator ICq. Following Remark 3.4 we know that 
/Co : L^(iM) -^ i^(jE) is of trace-class because it is the composition of HS operators. It should 
be also clear that the Fredholm determinants of Id-y^^ — K. and IdiR — /Cq are the same. Indeed, in 
the series that computes the two determinants, the contour 7^ can be continuously deformed to iR 
given the analyticity of the kernel. We now conjugate the integral operator on L^(zM) defined by 
/Co by the following Fourier-like unitary operator 

v2nTJtR V^irr Jr 

Consider now (the kernel of the operator defined by) To ICq o 7~~^: each of the terms in its defining 
sum in eq. (3.17) gives 

Tog,oToT-'=:Caix,y)^j-^ f ^e^(--) /d/x /dA ^-^^^^^e^^ (3.19) 

Now, if a; > a we can close the ^-integration with a big semi-circle in the right half plane, picking 
up only minus the residue at /i e 7^; viceversa, if x < a we close the ^-integration with a big semi- 
circle in the left half plane, which yields zero since there are no singularities within this contour of 
integration. In summary 

Iff C 3 xfj.+yX 

d/i / dA — X > a 



Ca{x,y)^{ (2^*)' Ah J^R (a^ " ^) =Ka{x,y) (3.20) 

X < a 

This proves that (using Remark 3.1) that 

ToJCqoT-^ ^Km\j (3.21) 

is the Airy kernel restricted to the union of interval /, concluding the proof. Q.E.D 

Looking at Theorem 2.1 we verify that in this case H{M) = where M(A) = 1 + G(A) as in Def. 
3.1; indeed the vectors f,g (3.5, 3.6) satisfy the stronger identity f^{X;d) ■ g{fj,;d') = when A,^ 
both belong to the same 77? (or 71,), respectively; indeed this implies that the additional integrands 
defining H{M) of formula (2.5) are identically zero. Thus we deduce immediately the following 

Theorem 3.2 The Fredholm determ,inant det(Id — Kxi), with I = [01,02] U . . . U [aAr_i,ajv] for 
even N or I — [ai, 02] U . . . U [on, 00) for odd N is equal to the isomonodromic tau function of the 
RH problem (3.3) namely 

da, lndet(Id - i^AiX/) = WM(5aJ:= / Tr(rii(A)r'_(A)Sa,,(A)) ^. . (3.22) 



10 



Theorem 3.2 implies also some more explicit differential identities by using the Miwa-Jimbo-Ueno 
residue formula; note first that the jump matrices M can be written as 

M(A;a)=e^(^)Moe-^(^) (3.23) 

where Mq is a constant matrix (consisting of only ±1 and 0) and 



T(A;a)=diag(ro,Ti,...,Tjv) (3.24) 

N 

To: 



1 ^ 

— ^i?,^ T,:=n-^ae (3.25) 



N . ^ 

The matrix ^(A; a) :— r(A; a) e^^^'""' solves a RHP with constant jumps and hence is (sectionally) 
a solution to a polynomial ODE. It was shown in [8] that (adapting to the situation at hand) 

Tr( -'^^ 



:{rz\xr_{\)EoJX))—=-rcs^Tr{r-\X)r'{X)da,T) (3.26) 

We then find 



Proposition 3.2 The Fredholm determinant det{ld — Kxi), with I = [oi, 02] U . . .U [ajv-i,aAr] /or 
even N or I — [ai, 02] U . . . U [ajv, 00) for odd N satisfies 

da, lndet(Id - Kxi) = -ri;,+i,,+i, (3.27) 

where Fi := liniA->.oo X (r(A) — 1) and T(X) is the solution of the RHP in Def. 3.1. 

Proof. This follows directly from formula (3.26) and from 

AT(A;a) = -A(^l-i..,,,,,) . (3.28) 

Indeed, plugging (3.28) into (3.26) we find 

da, lndet(Id - Kxi) = j^^ Tr Ti - Ti,t+ii+i- (3-29) 

Since detr(A) = 1 it easily follows that TrFi = 0, whence the proof. Q.E.D 

Now let's consider the case A^ = 1, i.e. the case in which we study the Airy kernel on the 
semi-infinite interval [s, 00). In this case it is immediate to observe that the RH problem 3.3 is 
nothing but the Riemann-Hilbert problem for the so-called Hasting Mc-Leod solution of the PII 
equation 

q"{s) = 2q^is) + sq{s) 

(3.30) 

q{s) ^ Ai(s), s — > -l-oo 

11 



(see for instance [18]; actually in the standard form the phase should be dxifJ.) = «-|- +ixfj, but the 
two different formulations are equivalent). Using our approach we can give a proof of the connection 
between the Airy kernel and Painleve II equation, as stated for the first time by Tracy and Widom 
[31]; note that this proof has purely academic relevance inasmuch as it is of different nature from 
the original one. 

Corollary 3.1 Consider the semi-interval [s, oo) = I, N = 1; then we have that 

\ 

{x- s)q{xfdx\ (3.31) 



det(ld - iCfs.oo)) = exp I - / 



where q{x) is the Hasting Mc-Leod solution (3.30) of the Painleve II equation. 
Proof. The RHP 3.3 for A^ = 1 implies the symmetry 

r(A) = air(-A)ai , (3.32) 

since the jump matrices have the same symmetry. In particular 

i^i(s) :== lim A(r(s,A) -1) =p(s)cr3 + ig(s)cr2 (3.33) 

A— >oo 

On the other hand we have already proven that the determinant is the isomonodromic tau function 
and from Proposition 3.2 we find 

a,logdet(ld-if[,^oo))=p(s). (3.34) 

Now it is well known and easy to prove (see again [18]) that the matrix VE'(A) :— r(A) exp (^z? 2,0-3)) 
satisfies a RH problem with constant jump and, consequently, the Lax system 

a^*(A;s) = C/(A;.s)*(A,a;) 9a*(A; s) = F(A; s)*(A; s) (3.35) 

with 

t/(A; s):= -^^3 - q{s)ai , V{\- s) := -A[/(A; s) + L' [s) - ^) ^^3 + V(s)^2 (3.36) 

Using the compatibility condition between the two equations above one can show that p'{x) = 
~q^{x) and q{x) solves the PII equation. Integrating twice we get the formula (3.31) while the 
asymptotic of q{x) can be deduced from its Stokes parameters as done in [22] and recalled just here 
below. Q.E.D 

3.1 Asymptotic behavior of the Hastings— McLeod auxiUary matrix 
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In the following we need some information on the behavior for a — > 
+00 of the AiO-RH problem (3.3) with phase ^„(C). The AiO-RH 
problem for N=l is nothing but the well known RH problem for the 
Hasting-Mc Leod solution of PII as discussed in Corollary 3.1. We 
need some simple estimates for the matrix Y A:i{C,\c!^ as ct — >■ +00 
and Id -^ 00. To this end we introduce the scaling ^ — y/az. The 
saddle point for the phase functions are: 



^9<x(^/CTz) ='^M ^ 




3 -^ zc = ±l d, = T\cri. (3.37) 

Figure 1: The jump matrices 

The contours of the jumps can be continuously deformed so for the Hastings-McLeod solu- 

, ii, 1, j-T, JJ1 -J- _Li -(• 4-u ^, tion of PII in the rescaled z- 

as to pass through the saddle points zi 2 := ±1 01 the phase on 

plane, 
the regular steepest descent contours 5 {z^ /i — 2:) = 0. It appears 

from (3.37) and the shape of the jumps (Fig. 1) that both jump 

matrices are uniformly close to the identity in any U" norm (including p = 00) with all these 



norms being 0(e 



by simple steepest-descent estimates. Since the jumps are analytic and 



the norm estimates hold also for deformed contours going to infinity within suitable open sectors 

(|arg(z)| e 



I + e, f — e) for 7/5, for example) the small norm theorem guarantees (for large a > 0) 

l|r.,K;.)-i||<^<£#S^. (3.38) 



1 



M 



1 + ICI 



with C some (inessential) constant, independent of ct, Q. What will matter for us is that the ct- 
dependence of T Ai is exponentially small as ct — > +00. 

Remark 3.6 It is important (but also simple) that the contours for the jumps of the Hastings- 
McLeod '^ -functions can be deformed to the rays of angles ±7r/4. 



4 The Pearcey kernel 

In this section we perform the same analysis as the one for the Airy kernel for the case of the 
Pearcey kernel 

eex(p)-e^(A) 



Kp{x,y,T) := 
e,(/x) := 



1 



(27ri)2 



l^^ T 2 



d^ / dA- 

7iU7R JiWL A — fl 



(4.1) 
(4.2) 
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The contours ^l,1r. are indicated in Fig. 2. Our aim is to 
consider the Fredholm determinant deiilA—KpXi) where we denote 
with the same symbol the Pearcey kernel and the related operator 
acting on L^(IR). xi is the characteristic function of the collection 
of intervals / := [ai, 02] U [03, 04] U . . . U [a2Ar_i, a2jv]- The relevant 
RH problem in this case is the following: 

Definition 4.1 Giuen /= [01,02] U. . .U[a27v-i7a2Af] we define the 
related PO-RH problem (Pearcey operator Riemann-Hilbert prob- 
lem) 

r+(A)=r_(A)(l-G(A)), Ae7LU7flU«, 




r{x)=i + o{x-^), A— ^00 



Figure 2: The contours for the 
Pearcey kernel and Riehrann- 
Hilbert problem (4.4) 



G(A): 



e'S'''i(^) -e^-2W 




(„)2Af+lge„2„(A) 





X7RU7, w+ 





e-e<.i(A) 



XimW 



(4.5) 

which consists of finding an analytic matrix-valued function T(X) on C\ (7^ U 7/j U iM) normalized 
to the identity at infinity and such that the limiting value ofT(X) approaching the contour from the 
left and from the right (r^{X) and r_(A) respectively) are related though the jump matrix G(A). 

Remark 4.1 The (— )2^-i (which -of course- is —1) is only reminding that the signs in the first 
row of G{X) are alternating. 

Proposition 4.1 The Riemann-Hilbert problem (4-5) is the RH problem associated to the integrable 
operator 

/^(A).g(M) 



K{X,fi) 



A — /i 



/(A) := 



1 

27ri 



/ 











X-ylUjrW 





-ieo(A)+aiA 



\ 



X«(A) 



5(A) 







gieo(A)-aiA 
(•_^2Ar+lgieo(A)-a2„A 



X-ilUjrW 



- jeo(A)+a2NA 

e" 5610(A) 





(4.6) 



(4.7) 



Xjr(A) 



(4.8) 
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Proof: It is just enough to verify that l-2z7r/(A)(f(A)'^ = l-G'(A). Q.E.D 

Theorem 4.1 The following identity holds 

det(Id -k)= det(Id -KpXi) 

where the operator defined by K is of trace-class on L^ijL U 7^ U iM.). 

Proof. The proof is essentiahy identical to that of Thm. 3.1 being, in fact, simpler since there 
shall be no need of contour deformations. Writing out K in full we have 

2i7r A — /i 

Then, using matrix notation subordinated to the split L'^ {jl^Ir)® L'^ (i^) we can write det(Id — A") 
as 



det 







J" 



EZi-y^'Sa, 



2N 



=.det[Id-^(-l)J"+ig„^. oj-] 



(4.10) 



where -this time- we have defined the operators J^, Qa^ as 



Ga -.L^ijL^-fR) 



fif^) 



e-5©o(C)+< 



f-> 



T:L^{iR) 
3(A) 



27rz 
i'(7fiU7L) 



'■/L'J-fR 



gieo(/.)-a/./(M)_du 



H^- 



e-s^of^^^^dA. 



(4.11) 



(4.12) 



27ri 7jR ^ /i - A 

It is a simple verification that the operators Qa,J- are of Hilbert-Schmidt class on "H :~ L'^{jl^Jr)® 
L^(ilR); in fact the reader may verify exactly as in Thm. 3.1 that both are of trace class by adding a 
vertical line to the left/right of the imaginary axis and writing them as composition of HS operators 
(the check needs to be done separately for 7^, 7/}). Note that the operator Id — X];=i(^l)"'^^^aj °J^ 
appearing in the last equation of (4.10) is an operator acting on L^(ilR) into itself and it is also 
of trace-class because resulting from the composition of HS-operators on H. We would like now 
to conjugate this operator by the (bounded) multiplication operator M := e^s®''^^^ on L^(«E): 
this step needs some clarification, since the inverse of this operator is unbounded. By inspection 
of the operator Qa we see that it automatically produces a result in the image of jM and hence we 
can multiply it on the left by A^^^. On the other hand, multiplying on the right by M causes no 
problems. The kernel of the resulting operator is thus: 

g©o(M)-0o(A)-a/i 



7W-iog„o J-oX(e,A) 



(27ri)2 



oCa 



d/i- 



71-U7R 



($-a*)(m-a) 



(4.13) 
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Note that one may verify directly that the kernel in (4.13) defines a trace-class operator on L^(iM) 
by realizing it as the composition of HS operators likewise as above. The resulting operator is 
obviously still of trace-class so that the Fredholm determinant is well defined (and has the same 
value as the original one). Now we consider the Fourier transform 

f{^i) ^ -L= f e^^/(^)d/x h{x) K^ -L= f e-^^/i(x)da;.^'^--^^) 

V2i7r Jm V 2i7r JiR 

Using (4.10) it is enough to prove that T o M~^ o {Ga2j-i - Qa2j) oT o MoT^^ = KpX[a2,^ua2j]- 
The kernel of the operator ToA^^^o^^oJ^oAlo T^^is equal to 



Ca{x,y) := —-, / -^e«('^-^) / d^ / dA — -, — e^^ (4.15) 



_L^ / ^e^''^-^) / dM / dA — -, —-y^ 

(27ri)^ Jai 27ri J^^vj^a JiR (^ - NiP - ^) 

Now, if X > a we can close the ^-integration with a big circle on the right, picking up only minus 
the residue at ^ e 7^; viceversa, if a; < a we pick up the residue at /z G 7l with the opposite sign; 
namely 

7- — -77 / du / dA — ; — — X > a 

£aix,y)=\ (4.16) 

7^^-^ / d^i / dA — — x<a 

. (27ri)^ J^^ JiR in - A) 

Formula (4.16) should be compared with formula (3.20). Thus we have 



1 /• f Qe.{t^)-ey{\) 

d^ / dA — — r — = Kp{x,y) x £ [02^-1, 02^] 



Ca2,^,{x,y) - Ca2,[x,v) = \ (2^0^A.u,Hy« (Am) " ' ' (417) 

I X ^ [a2j_i,a2j] 

and this last equation concludes the proof. Q.E.D 

From Theorem 2.5, since we have H{M) = by the same argument used after the proof of Thm. 
3.1, we deduce immediately the following 

Theorem 4.2 The Fredholm determinant det(Id — Kxi), with I = [01,02] U . . . U [a27v-i, a2Af] 
is equal to the isomonodromic tau function of the RH problem (4-5) (note that all intervals are 
bounded in the case of the Pearcey kernel). 

da,lndct{Id- Kxi) = iOM{dy.= [ Tr(rii(A)r'_(A)Sa(A));^ (4.18) 

J-rnUyr. ^ ' 271-1 

Moreover 

da, lndct(Id - K^xi) = -ri;£+i,^+i. (4.19) 
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Proof. The proof is identical to the case of the Airy kernel in Thm. 3.2 and Prop. 3.2, simply by 
changing the phase function Qa from that of Airy (3.1) to that of Pearcey (4.2). Q.E.D 

4.1 The solution for a = b 

When iV = 1 and a = ai = bi = b the Fredholm determinant is identically unity, since the action of 
the kernel is trivial. The Rieniann-Hilbert problem is -however- not immediately trivial. A closer 
inspection reveals that the solution can be written in closed form in terms of Cauchy integrals as 
follows: if we call A, B, C the three columns of F, the jump on 7 :— Jl^IRj i^ a-nd the asymptotics 
imply 



B{\) = 



A{X) 



1 

2i7r 

1 

2i7r 



s'3aMA{^i)d^l 



C{X) = 



II — X 



1 

2i7r 



'1' 




"0' 





+ 


1 







1 



and hence the solution is explicitly written for any a e C as 



r(A) 







.-e„ 



{fi — X)1m 



e®°(^)^(^)d/i 
/i — A 

/■ e-Q°(^)de 

JiR C ~ A 
d/i 



1 

2i7r 



.e„(M) 



(4.20) 



(4.21) 



(/I — A)2i7r 



{i - X)2m 

e-Q°(^)dg 
{i - A)2i7r 



,e„(A<) 



-e„(5) 



d? 



(fi — X)2m 



~ («-M)2i77 



f 



d/z 



d/i 



/ 



:.e„(/.) 



-ea(e) 



{fi — A)2i7r 



di 



^dM 



JiR (5-/i)2i7r 

(/i — A)2z7r 



d/i 



(4.22) 



/^ (/i-A)2i7r 

This is consistent with the triviality of the Fredholm determinant (identically one, since the kernel 
is null) , which suggests the solvability in "trivial form" of the associated RHP. Of course a similar 
expression holds for an arbitrary union A^ > 1 of "empty" intervals. 

4.2 A Lax system for the Pearcey Process and PDEs for the gap proba- 
bility 

Using the Riemann-Hilbert problem (4.5) it is possible to deduce a Lax system for the Pearcey 
process. We will treat the case of a single interval with endpoints a, b. Let's introduce the matrix 

^{X;a,b,T) :^ F(A; a, 6, T)e^(^'"'''^^) 



T{a, b, r; A) 



^diag(ea(r; A) + Qtir; A), -2ea(T; A) + efa(r; A), -2et(r; A) + Qair; A) 



ea(r;A) ^^-^X'-aX. 



(4.23) 
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In [2]^ two PDEs were found for the Fredholm determinant; we will show that these equation 
and a third one can be found in this framework, which therefore shows integrability by providing 
a Lax system. We will state the equations to be verified and explain the logic of the verification; 
the actual computation involves a significant amount of completely straightforward algebra and it 
is best handled by a machine. 



Proposition 4.2 The logarithm of the Fredholm determinant g{a,b,T) := log dot (Id- 
satisfies the differential equations in Oe '■— da + db, dr, e := ada + hdb 


'^^Ml) 


d%g + (S{dlgf-^Tdlg + 12dlg = Q 


(4.24) 


(-3e - 2Tdr + 2drdl + 1) dsg + l2{dlg){drdEg) = 


(4.25) 


e{l2drg ~ 2dlg) + (^d^g + ^d^dlg - Ad%g - 8 (9|g)') t+ 
+4a|5 + 16 {dlgf + 8 {dEd.g) d%g + 10 {dlgf + 16 {d%g) d%g+ 
+d%g - md^^g + 4dldlg - 24 (dEdrgf - 8 (drdlg) d^g - 85,3+ = 


(4.26) 


The latter equation (4.26) is new. 





Remark 4.2 By taking derivatives of the first two and adding them to the third with appropriate 
coefficients one obtains a somewhat shorter equation 

d%g - 8drg + 4r {2dlg - d%g - 2{dlgf) + I2e(drg)+ 

+l&{dlgf + Adldlg - 2A{drdEgf + I6d%gd%g + W{dlg)^ = (4.27) 

It is apparent that this equation cannot be obtained from the first two because it contains a first-order 
derivative w.r.t. r which is absent in the the other equations. 

Remark 4.3 If we had only the two equations (4-24, 4.25) the general solution would depend on a 
functional param,eter. The third equation makes the general solution depend on a finite number of 
parameters; this fact is a manifestation of the underlying isomonodromic system. 

We indicate how to verify these statements, but the actual verification (in particular of the 
third equation) was performed by a machine (the file is available upon request); first we use the 



"In order to obtain exactly the same equations it is necessary to rescale t i— >■ r/2 
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new variables 

i? = ^ , W:='^ , dE = d,+dt, dw^da^dt, e:^EdE + Wdw (4.28) 

and use the same symbols for the function ^ as a function of the new variables -E, W. The matrix 
"^{E, W, t; A) solves a RHP with constant jumps and has unit determinant; standard arguments al- 
low to conclude easily that (the dependence on the variables E, W, t shall be understood throughout 
for brevity) 

9^*(A) = i^(A)*(A) dE-i'iX) = i£;(A)*(A) dw'i'W = Lw(A)*(A) (4.29) 

*(A)'=^(A)*(A) (4.30) 

where the matrices Lt.LetLw and A are polynomials in A of degrees 2, 1, 1 and 3 respectively 
and / denotes the derivative in A. 

They can be found in the spirit of the [25] as follows. Denote by Vj the coefficient matrices in 
the (asymptotic) expansion of r(A) near infinity 

OO p 

r(A) = l + 5]^ (4.31) 

In this section we shall then understand F as the formal series appearing in (4.31) and each instance 
of r~^ shall be understood as the inversion of the above formal series in the sense of formal series. 
Since we know that all the matrices Le^Lw ,Ltt A are polynomials, they can be found by taking 
the positive part in powers of A (including the constant) of the following expressions 

Lr{X) =fdrTf-'' + drff-^ L e{\) =fdETf-^ + dETT^^ (4.32) 

Lw{\) = fdwTf-^ + dwff-^ A{\) = fT'f-^ + P'f-i (4.33) 

The important point to make is that the terms like Q^rr^^ etc. are formal series with only negative 
powers, hence, for example 

Lr = (fdrTf-^) , A = (tT'T-^) (4.34) 

where ()+ denotes the part with nonnegative powers of A. It should be evident that this is a 
polynomial A but also (and more importantly) a polynomial in the coefficient matrices Tj of the 
formal expansion, and involving only the first few (the first 2 for Lr and the first 3 for A). 

Consider now the ODE in the spectral variable A rewritten for the matrix F (and also for its 
formal expansion F) 

F' + Fr' = AF. (4.35) 
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This is at the same time a ODE that r(A) itself satisfies and also its formal expansion F: in terms 
of the coefficient matrices Tj this is an infinite system of polynomial equations in their entries. The 
coefficients of the positive powers of A in (4.35) satisfy identically the equation by the definition 
oi A = (rT'r^^j . The coefficients of the negative powers yield an infinite set of polynomial 
relations between the matrices Tj. 

Consider now the equation w.r.t. one of the parameters £', W, r (let's take the example of £'); 

dET + TdET = LET <^ dET = LeT-TOeT (4.36) 

A straightforward computation yields 

OeT^ ^diag(2, -1, -1) =: ^k, (4.37) 

From this we have L^ — ^xi and L^ — ^ [Fi,>fi] by looking at the coefficients of A^jA''. The 
coefficient A^-' yield 

dETj = ^ [^i,F,+i] + ^ [Fi, xi] F, , J > 1. (4.39) 

A similar computation for dr yields 

5rri , A2 / , Fi , F2 



i + ^ + ^ + ...Ui= 4'^a2 + lWa + 4") 



ri F2 
T + A^ + - 


. . V40) 




(4.41) 




(4.42) 



A 6 V A A2 

Lr^ - (-xiA^ - A [Fi, xi] + [xi, F2] - [Fi, xi] Fi) 

drTj^ ^ ([r,+2, xi] - [Fi, xi] F,+i + ([xi, F2] - [Fi, xi] Fi) F,) 

Equations (4.39, 4.42) allow to express any derivative (of any order by using Leibnitz rule) of the 
coefficient matrices T j in terms of polynomials in the same coefficient matrices (the expressions 
become larger and larger but stay finite at each step). In order to verify Prop. 4.2 (and eqs. 
(4.24,4.25)) it suffices to write 

OeQ^ - res Tr (F-^F'^i^r) = \ Tr(FiXi) , (4.43) 

A— 00 o 

5,g= - res Tr (F-IF'^.T) - ^ Tr ((F? - 2F2)xi) , (4.44) 

A— 00 O 

One then makes repeated use of eqs. (4.39, 4.42) plugging the result into equations (4.24,4.25); the 
result is a rather large polynomial in the coefficients matrices Fi,F2,F3. 
One then compares with the following polynomials 

res Tr((F'F-i+FT'F-i-yl)a£;T)^ , i = l,2,3 (4.45) 

A— 00 A-^ 

which have to be identically zero due to (4.35) and verifies that they are the same polynomial thus 

completing the check. 

20 



5 Large asymptotics for the Pearcey Gap probability and 
factorized Tracy— Widom asymptotic behavior 

The goal of this section is to prove the asymptotic factorization of the Fredholm determinant for 
a "large" gap in the Pearcey process into two Fredholm determinants for semi-infinite gaps of the 
Airy process. In detail 



Theorem 5.1 Let Kp(x,y;T) denote the Pearcey kernel and KAi{x,y) the Airy kernel. Let 

a_ =a_(p) := -2A^ + Apv^ a+ = a+(cr) := 2A^ - Act v^ r := 3A^ (5.1) 

then OS A — )■ +00 



det Id-i^p(»,»;T) 



det Id- KAi{»,») 



[a^M+] 



det Id-/fAi(»,») 



[p,oo) 



(l + 0(A-i)) 



[cr.oo) . 



and the convergence is uniform over compact sets of the variables p, a of the form 

K. 



Ki<p,a< ^A^ -00 < Ki, K2 < 2, 

73 



(5.2) 



namely for which p, a are uniformly hounded below (w.r.t. K) and do not grow faster than r 3 . 

More generally when we have a union of intervals; 
Theorem 5.2 Let 3 be a union of intervals of the form 



(5.3) 



J K 

'3^ Ij[a2i-i,a2j] U [a2j+i,6o] U |J [62/0-1, &2fe] 
]=i k=i 

ae:= a^ipe) = -2A^ + Api^^ be -.^ a+{cr2K+i-e) ^ 2A^ - ^cr2K+i-e\/i , r = 3A^ (5.4) 
Then, as A — > cx) we have 



det [U-K^{;»;t] 



det Id-XAi(»,») 



det Id-XAi(»,») 



Ji, 



(l + 0(A-i)) (5.5) 



J2, 



K 



Ji — U [^2.-1 ' P^t] U [P2J+1 , 00) J2 := U K,_i , cr,J U [a^^^^ , cx)) (5.6) 



t=i 



In the above it is allowed that a2j+i = &o (hence the middle interval is missing) in which case 
Ji ■=Ue=iiP2i-i,P2i] and J2 := UfciK^-i, cr2J- 
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t = 1 



Figure 3: The Pearcey and Airy processes. 

The proof of this theorem consists of the whole present section; it rehes essentially upon the 
Riemann-Hilbert problem constructed in Sect. 4 and the Deift-Zhou [16] steepest descent method. 
We will start with considerations that apply to the more general case of Theorem 5.2 but then 
specialize to the case of Theorem 5.1 in order to avoid unnecessary complications (which are purely 
notational and not conceptual). 

Remark 5.1 The parametrization of the endpoints a^,a-\- in Thm. 5.1 (and of aj,bj in Thm. 
5.2) has the following meaning. The Pearcey process arises in the study o/ self-avoiding random 
walks on the line, conditioned to start at the same point (say the origin) at time t = Q and end at 
time t — ti > at two distinct points moving away from the origin as N^''^ (N being the number 
of particles); at any time < t < ti the bulk of the walkers consists of either one or two finite 
intervals. There is a critical time < tc < ti at which this bulk undergoes a transition from a 
connected interval to two intervals. The two new emerging endpoints [a(t),b(t)] move away from a 
common point a{tc) — b(tc) according to a{t) — ac + 0{t — tc)^ . The Pearcey point process describes 
the statistics of the random walkers in a scaling neighborhood oft — tc and a,b = Oc', more precisely 
we rescale t >->■ tc + N^^/'^t, a ^^ Oc + N~^/'^a and b t-^ Oc + N^^/^b; see for instance [33] . The 
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asymptotics as t = 3A^ — > cxd and a — a-,b = a^ as given in Thm. 5.1, is the regime where we 
look "away" from the critical point and it is expected to reduce to two Airy point processes, which 
describe the edge-behavior of the random walkers. 

The phase 6a (t; A) (4.2) has an inflection point with zero derivative when the discriminant of 
the derivative vanishes: 

DiscrimA(A''' - rA - a) = -At^ + 21 a^ = (5.7) 

The neighborhood of the discriminant is parametrizable as follows: 

a± (s) := ±(2A9 - As^) , r = SA^ (5.8) 

We have the expression 



e.,(r; A)= ^A^^ - sUSt (^C^ - <) +i^±^ = (5.9) 

= C(s)Tt?.(C)+^^^^ (5.10) 

C:=CL,fl = 3U(A±A3) (5.11) 

For definiteness we consider only the case a2j+i < &o- 

Preliminary step: we conjugate F by the diagonal constant matrix Dq with entries 



, /2,7+l 2A'+1 \ 

— T Y.^M+ E ^M (5.12) 



{Do)ii ■■= 2^ 

{Do)e+i,e+i := iDo)n - Cipe) , l^l...2K + l (5.13) 

{Do)i+2.j+2,i+2j+2 = {Do)ii-C{(Jt) , l=l...2J + l. (5.14) 

with C{s) = |A^^ — 33 A^s as in (5.10). We denote by a tilde the new matrix and respective jumps 

r(A) := c-^°F(A)e^° , M(A) := e--°"M(A)e^° . (5.15) 

The conjugation by Dq has simply the effect of replacing the phases Qa^{pi)i 0a_((Tj.) by ^a^ipt) — 
C{pi) and Oa+((Tj) ~ C{(Ti), respectively, so that their critical value is zero. 

Factorization of the jumps: the jump on the imaginary axis Mq := M\ for F can be 
factorized 

Mo = Mo,rM-1 , Qa^is)W ■■= e,±(,)(A) - C{s) (5.16) 
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10,B 
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Figure 4: The contours of jumps for the RHP for Y. 



Mo,. 
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l2,7+l 








'i-2K+l 



Mo,_ 



1 






-e"®"-'"!' 


12J+1 




-C^^^-Csj+l) 










12a:+i 


L 







(5.17) 



This factorization aUows us eventually to replace the jump on iR by a jump on two separate contours 
in the right/left halfplanes; this will be used in the final step below. 

Rescaling: in view of the fact that the critical points A = iA"^ escape to infinity, we now 
introduce the scaled variable z := A~^A. 

Final step: We shall choose the contours for our final matrix as follows in the plane z = A^'^A: 

• Two counterclockwise circles centered at z = ±1 of radii r — A^^ (i.e. of radius 1 in the 
A-plane) ; 

• the two contours 'Jr.l shall be moved to the contours consisting of the two pairs of straight 
half-lines originating a z — ±1 (A = ±A^) with slopes ±^; 

• the two contours that shall support the jump that was originally on zM will be denoted by 
70, fl and 70, l = -7o,fl- The contour 70,7? is 



7o,fl := {3?z = 1 - A^^e^ , |z - 1| > 1} U {1 + se=^T^ , s € [0, A'^)}, 



(5.18) 



24 



and oriented upwards. 
The arrangement of these contours is best illustrated in Fig. 4. Define now 

{r{A^z)Mo^R{A'^z) in the region between zM and 70, _r 
r(A^z)Mo,L(A^2:) in the region between zK and 7o,_r 
r(A''z) elsewhere. 

The matrix Y{z) solves a new RHP {M{z) = M(A^z)) with jumps: 

Y+^Y^M , Y{z) = l + 0{z-'^) , 
M{z):^ Mo^ii{z)xjo_a + Mo,l{z)x^o.l + MR{z)xjn + Ml{z)xil 



(5.19) 



(5.20) 



Ml = Mr := 



f 


_C§a_(Pl) ^ ^ _ (_)2J+lg§»-(P2j+l) 


e®"+(-2if+i) . . . (-_-)2Arg§a+(^j) 




I2J+I 








l2i<:+i 



(5.21) 



Mo._ 



L2,7+l 



"+^"2^ + 1' 



e-e^+c-i) 



'-2K+1 



Ml 



0,R 



1 






-e"^"- ("2,7+1 > 


I2J+1 









I2K+1 



(5.22) 



e,^(,)(z) := e,^(,)(3A6; A^z) - C(s) = T^,(C) + A^^ ^^^ 



(5.23) 
(5.24) 



(The subscripts L,Rhi (^ mean that one is centered around the Left critical point z = — 1 and the 
other around the Right critical point z = 1). The two circles for the time being do not support any 
jump of Y , but will be used eventually in the construction of the approximation. 
The main idea of the proof is now that 

• The matrix Mq.a is exponentially close to the identity in L^ijo.R) p < 00 outside of the disk 
at 2 = 1. Similarly for Mq.l relative to 70, l ■ 

• For the matrix Mr, the entries of the form 6"°+'"*' are all uniformly small (as A — >■ 00) in 
any ^^(7^^), p < 00; ditto for Ml and the entries e""-'"*^'. 
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• Inside the disks we will use an parametrix build out of the AiORHP (3.1) of the appropriate 
size that solves the jumps up to exponentially and uniformly small terms. 

The proofs of the first two bullet points rely upon Lemma 5.1, Lemma 5.2 and Lemma 5.3. 

Lemma 5.1 Let K2 < 2 be fixed and a < K23^^A^; then the function e"°+<"''^' (t — 3A^j tends 
to zero exponentially (in Aj in any norm, LP(-fji), 1 < p < 00 

lle^"+<-'^^)|U.(,..)<e-^(^-^^)^". (5.25) 

Similarly, the function e""-''"'^^-^ tends to zero exponentially (in Aj in any norm L^^^l), 1 < p < 00. 

Proof. We set z = 1 + (1 ± i)t, t E M+ as a parametrization of 77^. A trivial computation yields 
(for both signs!) 

Hit) := 5R(e,^(,)(z)) = Ai2 l-t^ _ 2^3 + J ^ - 4 j j + 2A4 (a^ - 2A8) (5.26) 

We have |e""+'"'| = e^^*-': under the assumption that a < K2S^^A^ with a constant K2 < 2 
(independent of A) we have (recall that t e M+) 

iJ(i) <-Ai2t4-2Ai2(2-i^2) (5.27) 



|e^|U.(,,) < e-2(2-^^)^^^ ( 2 I e-P^''''^dt ] ^ < e-'^'-^^^^'' ( ^ ..'L.^,.. ) ' (5.28) 



y2piA3r(3/4)/ 



One verifies that this last expression of p is less than 1 (for A > 1). Note that if K2 = 2 and a 
saturates the upper bound, one can still show that the L^ norms are exponentially small, except 
the L°° which is constant. The other case is completely analogous. Q.E.D 

Remark 5.2 In Lemma 5.1 the upper bound K2 = 2,(t = if23~^'^A^ corresponds to a_|_(cr) = 0, 
so the above estimate holds as long as a+ > and t — > 00 (see Fig. 3). 

Lemma 5.2 Let Ki > — cx) be arbitrarily fixed; then, uniformly in p > Ki the function e~°-(f'^^' 
tends to zero exponentially (in K) in any norm LP{'jii,e), 1 < p < 00, where "fB.,e '-^ IR \ {\z — 1| < 
A~'^} and specifically 

lle^-'^'^^^IUn^H,.) < e-^^^+^^l^^l^ (5.29) 

Similarly for the function e~°+<'''*-^'^ in L^i^L^e), lL,e '■= IL \ {\z + 1| < A^'^}. 



26 



Proof. The rays of ^n.e can be parametrized as We set 2; = 1 + e*'*' (i + A ''), t > 0. A trivial 
computation yields (for both signs!) 



,2 , .4 (^/2A3 + 2) t3 3 {V2K^ + l) ^ 
-(p)V )) I 4 2A3 2 A6 



//(i) := 5R(e,_(„)(z)) = A^M - f4 _ vv 7 _ IV ^^+ (5,30) 



/2p^ 2 + 3\/2A3\ \/2pv/3 l + 2\/2A3 



' \^ 2A« 2A9 / 2A" 4Ai2 ) ^^'^^^ 

The coefficient of the linear term in t is negative for p > K\ and sufficiently large A hence we can 
estimate 

H{i) < -A^^^ - ^A^ + ^A\K,m. (5.32) 

This shows that the L°° norm tends to zero since H{t) < iJ(0) < -^A^ + ^A|ii'i|v/3 ^ -00. 
As for the other L^ norms, we have 



e' 



^||P, -/ eP«(*)di<eP«(0) / eP(^W-^(°»di = eP^(o) / e"'^* di = e^^^ 



00 



7rV2 



Hence (for sufhcicntly large A the last bracket is smaller than one). 



|e^|U.(^,,j < e^(0) < e-^A^+#A|^il^ ^ Q (5.33) 



Q.E.D 



Lemma 5.3 Let Ki > —00 &e arbitrarily fixed; then, uniformly in p > Ki the function c ""-(p)^^-' 
tends to zero exponentially (in K) in any U" , 1 < p < 00, of the vertical rays of jo,ii, denoted by 
lo,R,e, in particular 

l|e-^"-<'"^^^llL.(,„,«,.)<e-^-^+^^^/2+i , 7oAe:=7o,i.\{k-l|<A-n (5.34) 

Similarly for e ~ <^+'.'='>^'^' ' on the vertical rays 7o,L.e -^ 7o,L \ {k + 1| < A^'^}. 

±2i7T/3 

Proof. Once more this is a trivial check using z = 1 + ^^^ — zLit, t (^ M+ 

i7(f):=-5R(e„_(^)(z))= (5.35) 

12 / 1 .4 1 /T^.ri / 3 3 \ ,5> 3v3,\ 1 1 



^ > 4 



< - V^' - A' - ^Ap + ^ < - V^^ - A' - KiA/2 + 1 (5.37) 



4 2^8-4 
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Thus we have 



)(-)| 



\LP{fo.R,±) 



= 2 



eP"('^dt 



< e 



_A3_^A+1 






2piA3r(3/4)/ 



(5.38) 



Hence it is clear that e-^^*^ satisfies the assertion (for large A, e.g. A > 4) because the constant in 
the p-th root is less than one for A large. Completely analogous computation on 70, l- Q.E.D 

In order to obtain the desired asymptotic information one needs to solve the RHP within the 
two disks centered at z = ±1; these local RHP are expressed in terms of (r.l (5.11); the local 
parametrices are then constructed in terms of the RHPs for the Airy kernel in the form (3.3). 

As we have anticipated, we will enter into details now for the simplest case where K = J = 
and hence we have a single interval for the Pearcey kernel of the form [a_(p), a_|-(cr)]. The jumps 
for Y are explicitly 



r(A)+-y_(A)Mo,fl 
^ " 1 



A e 70, fl, ; 
" 
1 



y(A)+ = r_(A)Afo,L Ae7o,L 

1 

1 



Mo,R := 



r(A)+-y_(A)ML Ae7L r(A)+-y_(A)Mfl Xe-tR 

Mr^Ml:= 1 

1 

Afl, :=(A-A3)=A3(z-1), 



Al:=(A + A3)=A3(z + 1) 



(5.39) 

(5.40) 
(5.41) 
(5.42) 
(5.43) 



Note that the matrices Mfi^M^ have identical form, but are defined respectively on ^r,^l- Also 
important is the fact that the terms A|,A|j (with AljA/j defined in (5.43)) appearing in the ex- 
ponents are uniformly bounded within the disks because these have radius A~'^ in the z-plane (or 
radius 1 in the A-plane). As a result of Lemmas 5.1, 5.2, 5.3, the jumps on the contours outside 
the two discs at itA^ are exponentially close to the identity matrix in any L^ norm (including L°^). 
Moreover, thanks to Lemma 5.1, the jumps on 7^^ and 7^ can be written as 



{z)= {1 + (c-2(2-^2)A^^)) M*')(^) ^ M 



Mn(z) = 



(0) ._ 

R -^ 



1 --eBp(Ci<)+\>-H 
1 

1 



Mr, 



M 



(0) 



.(z)=(l + o(e-2(2-^^)A"))M["), 
where the estimate holds uniformly on p, cr < 1^23^3 A® for any K2 < 2 



1 e"''-('^^)+J^^ 
1 

1 



(5.44) 



(5.45) 
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1 ;e-''^(a)+iAi 




MCl)- 



Figure 5: The exact jumps of Fai near 
z = 1. 



Figure 6: The exact jumps of Fai near 
z = -l. 



5.1 Approximation of the solution 

Notation: Given a 2 x 2 matrix A we will define the 3x3 matrix A^"^'^', 1 < i 7^ j < 3 as the 
matrix whose elements {i,i),{i,j),{j,i) and (j, j) are the elements (1, 1), (1, 2), (2, 1), (2, 2) ofyland 
zero elsewhere. 

5.1.1 Parametrix near z ~ 1 

In this subsection and in the following we use the Hasting-McLeod matrix FAi (defined in Section 
3.1 and whose asymptotics is discussed ibidem, and which is nothing but the solution of the RHP 
in Def. 3.1 for A^ = 1) as parametrix for the RH problem related to Y inside the disks centered at 
z = ±1. The same parametrix was used in several papers in different contexts, e.g. [6, 9]. Let's 
start with the disk centered in 2; = 1. Consider the matrix 



FAi(z) := [eS^«-^FAi(Cfl;p)e-^^«"^ 



(5.46) 

(5.47) 



It solves the jumps within the disk centered at A'^ as indicated in Fig. 5 (see Remark 3.6). Its 
behavior on the boundary A'^jz — 1| = \\r\ = 1 is 



FAi(^) 



AA^sy3 



pip) 



ie-* «(?(p) 



4A 



«9(p) 



-pip) 



+ C'(A-2) 



With that in mind we shall define the parametrix as follows 
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Vr{z) 

|z-l|=A-3 



1 



1 



-»(e"^ «-l)g(p) 

p{p) iq{p) 

-iq{p) ~p{p) 









TAiiz) 



0(1,2) I ^e ,, 



n (1,2) 
(5.48) 

(5.49) 



The notation ©(^'^^(A^^e^sP^ ) denotes a matrix of the indicated order only in the entries i,j £ 
{1, 2} and zero elsewhere (the bound as a function of p follows from the estimate (3.38)). We thus 
have 

Proposition 5.1 The matrix Vr fulfills the requirements: 

1. is hounded in the disk \z ^ 1\ < A^^ (\'^r\ < Ij/ 

2. within the same disk it solves exactly the jumps conditions with Mr replaced by Mj^ as defined 
in (5.44); 

3. on the boundary \\r\ — 1 has behavior (5.49); 

4. the determinant detT'i^ tends to one uniformly within the disk, and hence V^ is bounded in 
the same disk. 

Proof. Only the first and last point require comment while the remaining ones follow from direct 
manipulations of the definitions. The boundedness follows from the boundedness of Tai and the 
boundedness of the prefactor in (5.48); indeed the apparent pole for z = 1 4^ Xr = is not present 
since ^(e*i'^« — 1) are analytic at Xr — (i.e. z — 1). As for the determinant 



detVLiz) = 1 + 



1 2smh{X%/8)q^{p) 



A2 



3A|, 



(5.50) 



The function w^^ sinh(u'*/8) is bounded in the disk \w\ < 1 (which is our disk in the z-plane since 
Xr = A*(z — 1) and |z — 1| < A^"^). On the other hand the function q{p) is uniformly bounded in 
the range Ki < p < ^A* because q{p) is continuous on M (it is the Hastings-McLeod solution of 
PII) and decays to zero at p ^- +00 as the Airy function Ai(p). Thus the assertion (since there is 
a division by A^). Q.E.D 
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5.1.2 Parametrix near z = 1 

Similarly to the case z — 1, consider the matrix 



-iAl<T3 



rAi(z) := ^e5^^-^'a2-^rAi(CL;'T)a2e- 

Al:=A + A3 = A3(z + 1) , C^:=Al^A 

^0 -i^ 



0-2 



i 



(5.51) 



It solves the RHP with the jumps indicated in Fig. 6. Its behavior on the boundary |z+l| = A ■^ 
(lALl-l)is 



-p{u) — ie" ^(?(cr) 



_ 1 \4 



rAi(A) = 1 + xaTW 

With that in mind we shall define the parametrix as follows 





+ 0(A-2) . 



7'L(z):=diag(0,l,0) + 



Vl{z) 



1 



AA,^ 



4A 



-i{c J i - l)q{a) 



i[e~^^i. -l)qi<j) 




rAi(^) 



(1.3) 



-p{a) -iq{(T) 



iq{a) p{a) 



0(1,3),_, 



(5.52) 



|.+l|=A-3 A4(z + 1)^ 

The next proposition is proved as in the previous section 
Proposition 5.2 The matrix Vl fulfills the requirements: 

1. is bounded in the disk \z + 1\ < A^^ (\^l\ < Ij/ 

2. within the same disk it solves exactly the jumps conditions with Ml replaced by Mj^ (5.45); 

3. on the boundary \z + 1\ — A^^ (\^l\ — IJ has behavior (5.52); 

4. the determinant detVL tends to one uniformly within the disk, and hence Vj^ is bounded in 
the same disk. 

5.2 Approximation and error term for the matrix Y 

We will define 

<i>oiz) := 1 



F.(.)(-) , (.-F,(.)..)(M) _3<,,^, 



$(z):=<' "' ' ^A4(z-1) ^A4(z + 1) 

Vr^lW |^ti|<a-3 

where fi(s) was introduced in eq. (3.33). 



(5.53) 
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Lemma 5.4 For Ki < p,a < K23~^A.^ , Ki > — cx), K2 < 2 we have that the error matrix 
£{z) :— y(z)$~^(z) solves a Riemann-Hilhert problem with jumps on the contours indicated in 
Fig. 7 and of the following orders 



£+{z)=£^(z)Jr,l{z) 
Jl{z) = l+O'^^M 



\z±\\ 



A- 



(5.54) 



ae 3 
A2 



3 

0-2 



(1,2) 



Jr{z) = i+o^'-'H 



A4 



\+o 



U>5{1,2) 




whereas on the remainder of the contours we have 



£+{z) = £^{z)Jo{z) , Jo{z) = 1 + 0(A-°°). 



(5.57) 



and the big-O symbols are in any L^ norms (p < co). The subscript on Jq stands for "Other". 

Proof. On the boundaries of the disks the estimates (5.56, 5.55) (including the shape) foUows from 
(3.38) and (5.49, 5.52), the definition of $0 (5.53) and matrix multiphcation. 

The other jumps are on 70, _R,e, 70, L,e,7L,7fi; on 7o,fl,e,7o,L,e the 
jumps are ^oMq^r^q^ , ^qMqx^o^ ^^^ Lemma 5.3 together with 
the fact that both $07 ^0^ ^^^ uniformly bounded implies that the 
jump is within 0{e^ ^t^) of the identity in any L^, p < cxd. On 
7ij and outside the right disk the jump is ^qMr^q and now Lem- 
mas 5.1 and 5.2 imply similarly that the jumps are within 0{A^°°) 
of the identity. On 7^ inside the disks, since the jump of Vr is Mj^ 
(Prop. 5.1, point 2), the jump of f is Vr-Mr(M^^^)-^V 



(°'^-i-"j,l. Com- 
1 





Figure 7: The jumps of the 
error term £ . 



paring with (5.44) and using the boundedness of Vr^V]^ (Prop. 
5.1, point 4) we see that the jump is also exponentially close to the 
identity. Similarly for the jump on 7^ outside and within the left disk. Q.E.D 

The small-norm theorem (which we sketch in App. C) says that -uniformly on closed sets not 
containing the contours of the jumps which we collectively denote with 7- 

1 



ll'^W-'ll^diSi^'-'*-'* 



(5.58) 



where the constant £( J) is bounded by some multiple of the L^ norm and the square of the L^ 
norm of J{z) — 1 (here J{z) is the jump on all possible contours for £). Both of these two quantities 
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on a circle (or any compact contour for that matter) are bounded by a multiple of the length of the 
curve. Overall the main contribution to £( J) is coming from the two circles and since these shrink 
as A~'^ (so does their length) we have 

g:(J) < Const k^^\\J{z) - 1||l~ (5.59) 

with the L°° norm here above restricted to the two shrinking circles. The error matrix £{z) is then 
found as the solution of the singular-integral equation 

£(.)^l+/^-M(4H4)^ (5.60) 

J 2iTT{w — Z) 

and can be obtained by iterations 

J 2i7r(w — z) 

where the integral extends to all the contours of the jumps; however the only contributions that are 
significant come from the boundaries of the two disks (the other contributions give exponentially 
small terms). This way one promptly verifies that the matrix £{z) has orders in the following shape 
which follows from the particular shape of the jumps Jl, Jr in (5.55, 5.56) 



1 

dist(z. 7) 



g(^)=l+ ......,.. gp.^ 



O,,, := 0(1.2) /^g-|pfA _^ 0(1,3) /y|^-|.fA ^^fVPZ^-lial+ps) I (5g2) 



5.3 Conclusion of the proof of Theorem 5.1 

By applying Thm. 4.2 we have for a_ = a-(p), a+ = a+(o') with a± given in (5.8 



9a_lndet Id-i^^ 



-ri;2,2 , 9a+lndet Id-X, 



-ri;3.3 (5.63) 



where Fi := limA-j.00 •^(r(A) — 1); to compute the limit we take A within the sector j + e < arg(A) < 
■| — e (form some e > 0) so that the distance between z = XA^^ and the jumps of £ grows like 
0{z). then our approximation for r(A) from (5.62) 

r(A)=e^o£('AUJ^V-i^o 



A3 / " V A3 



\ ^ 0,,. ^ F^{p)i^■^) ^ {a^'F,{a)a,)(^'^) 



I ^ ^3A(A-A3) ^A(A + A3) ^ ' ^ ^ 
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Thus, computing the Hmit that defines the matrix Fi = HniA-^oo •^(r(A) — 1) we have 



Ti = e^° 0(i'2) 



\/pe~ 



A2 



From. (5.65) we have 



0(1,3) 



A2 



O 



/pac 



\(p3+a3) 



A2 



Fi(p)(1^2) (^-lj^^(^y^)(l,3) 



^A 



^iK 



^--^*'+°(^" 



O 



pcre 



^V2 



r.«^^PM + o 



A2 



o 



pcre 



\(pl+al) 



A2 



(5.65) 



(5.66) 



(5.67) 



Recall that p(cr) is the logarithmic derivative of the gap probability for the Airy process; noticing 
that da_ = 3^- dp, da^ = — ishirda, the above equations say (in the notation of differential forms) 



dp,<jlndet \ld- K^ 



d„lndet Id-^Ai 



Ai 



la-(p),a+{cr)], 



dp In det Id - KAi 



Ip.oo)^ 



+0(^e-t''Mdp + Ol— c 



dcr + O 






(dp + da) 



(5.68) 



where we recall that the estimates hold uniformly for Ki < p,a < ii'23"^^'^A^ and < K2 < 
2, Ki > —00. To obtain information in integral form we integrate (5.68) from some point (pojCTq) 
to the point {p, a) = A^/2(l, 1); this point is still within the range of validity of all of our estimates 
in Lemmas 5.1, 5.2, 5.3 (in particular Lemma 5.2) and corresponds to an interval [a_,a+] of fixed 
size (not growing as r = 3A^ -^ cx)) containing the origin in its interior. The correction terms in 
the differential form (5.68) yield a form that is integrable on any rectangular domain of the form 



i^i < P, 0- < — i:A^ -oo<Ki, K2<2. 
V3 

and the result is a function of order 0(A^^) thanks to (5.68). Thus we conclude 



(5.69) 



Indet Id - K^ 



j =lndet [ Id - i^Ai )+ Indet [ 



Id - KAi 



\+0{A-^) + C 



[p,oo) 



It is shown in App. B that the Fredholm determinant of the Pearcey kernel in this regime (i.e. 
when the interval [a_, a+] is fixed) tends exponentially fast to unity as 3A^ = r — > 00. On the other 
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hand also the Fredholm determinants of the Airy-kernels on the right-hand side tend to unity, and 
hence the constant of integration C must be zero. The proof of Thm. 5.1 is now complete. Q.E.D 

The proof of Thm. 5.2 is completely analogous, with the complication that the local paramet- 
rices would have to be built out of bigger-size version of the Hastings-McLeod auxiliary matrix 
constructed using Def. 3.1 and the analogue of the simple estimates in Section. 3.1 for this larger- 
size version would have to be derived. 

A Proof of Theorem 2.1 

We briefly recall the situation we already described in the introduction. Following [20] (and ref- 
erences therein), we start with an integrable kernel K{X,^) acting on L'^{C) where C is a analytic 
contour in C, possibly extending to infinity. Being integrable means that we can write the kernel 
as 

A — /i 
with /^(A)g'(A) = (i.e. the kernel is non-singular on the diagonal). Moreover we assume that 
vectors / and g are analytic in a neighborhood of C. Then the resolvent R = (Id — K)^^K is again 
an integrable operator related to the vectors F, G given by 

F{\) = r(A)/(A) , G(A) = (r-i)T(A)g(A) (A.l) 

where T is the solution of the RH problem (1.5)-(1.6)-(1.7). From now on we denote with a prime 
the derivative w.r.t. the spectral parameter. 
Using formula (1.4) we have 

aiog(det(Id -K)) = - Tr{dK + RdK). 

The first piece is immediately computable since ^(A, A) = f'^{\)g{X) so that 

- T^{dK) ^- j (a/'T(A).g(A) + f^{\)dg{\))d\ - H,{M). (A.2) 

For the other part we get, using the explicit expression of the resolvent, 

- Tr{RdK) ^1 l^^^-L^Tr [{fT-'){^,){rf){X}{df{\)f{fi)+g^{X)df{f,)))dXdf, 

= JIjx~^^' ((/5''r-i)(^)(r/(a.9T))(A) + {i^fWr-'){^^){rff)i\))d\d^^ 
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where we just used the invariance of the trace for cycHc permutations and, in the last passage, we 
used the two identities 



(r+-r_)(A) = -27r*(r/gT)(A) 



(A.3) 
(A.4) 



coming directly from the related RH problem. 

Now we work with the assumption that f,g are analytic in a neighborhood of the contour and 
we split the integral above in two part. The two parts will be singular when A = /i so that we have 
to pay attention to the order of integration. More precisely we write it as 



TriRdK) 



^ Tr ((rii - r;i)(M)(r/(95^))(A) 



c i>^~ f^) 
1 



dX 

27ri 



d/i 



dX 

2Tri 



dfi (A.5) 



We would like to apply Cauchy Residue's theorem but, in the first term, we have to interchange 
the order of integration. Given an analytical function F{X, ^) in a neighborhood of the contour we 
observe that we have 



_J: pa .dX 

c [X — /i)^ zm 



d/i 



'ciX^W^^^'^^^^ 



dX 



A=A4 



c(A^^^"''^ 



dX 

27ri 



d^j, 



M=M+ 



where in the last equation we just renamed the dummy variables. Hence using Cauchy's theorem 
and the equation above for the first term of the r.h.s. of (A.5) we get 



Tr{RdK) 



Tr ((r;irV/(9g)T)(^))d^- f Tr [{rz'T'_idM)ifi))di, (A.6) 



Hence, combining (A. 2) and (A.6), we get 



aiog(det(Id -K))^- TiidK + RdK) = 

H,{M)- I Tr ({T-'T'J{dgf + T-_'T'_{df)f){i,))diM 
J c 



(A.7) 



On the other hand we have that 



WM(a) = I Tr {v-jT'_dMM'^) ^^^ - J Tr (riir'_((a/)5'^ + f{dg^)){l+27rzfY)]dX 
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= -| Tr (rz'r'_{idf)f + fidf))+2mrz'r'_f{dgffY))dx 

= -j Tr{T-jT'_{{dJ)f + mf) + (r;i - rz'r_ f{dgr))dX 
= - 1 Tr (rz'rUdfW + ri'T'_f{dgf))dX (A. 



ic 
(here again we used (A. 4)). Now we observe that, because of the RH problem, we have 

r'_=T'+M-'+2mr+{ffy (A.9) 

so that 

y Tr(r;ir'_/>g)T)dA = /Tr hl'(T'+{l+27rtff)+27Ttr+{ffY)f{d9)AdX 

= J Tr (ri'T'+f{dgfyX-27Tz J TT{fff'{dg)^)dX (A.IO) 

where we used several times the invariance of the trace under cyclic permutations and the fact that 
f^{X)g{X) = 0. Substituting (A.IO) into (A.8) we get 

dX 



LOM{d) = /" Tr (TZ^T'_dMKr^^ 

=- J Tr (ri^r {df)g^ + T^'r'Jidgf) dX + 2m J Tr{ff f{dgf)dX = 

=- I Tr {TZ'T'_{df)f + T-+^T'J{dgf)dX + H^iM) (A.ll) 

J c 

Comparing (A.ll) with (A.7) we obtain (2.5). Q.E.D 

Remark A.l It can also be proven that 

H2{M) = — / TY{M'M~^dMM-^)dX 
2 Jc 

through the following computation 

f TT{M'M-\dM)M-')^= f Tr [{ffy{l+27Tiff)d{ff){-i-+27Tiff))dX 

= / Tr {ff{dfW+frmf)dX 

^ JTr {ffg^idfyfff'{dgf)dX = - / Tr {ffidgff + fff'idgf)dX 

= -2jjr[ffr{dgf)dX 

This same quantity is introduced in [8] and is nothing but the difference between um and the so-called 
modified Malgrange's form (see the cited paper for details). 
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B Behavior of K„ for r -^ +00 



We consider the behavior of the Pearcey kernel when t — > +00 and x, y belong to some fixed 
bounded interval [a,b]. The analysis is a regular case of steepest descent. 
We perform the rescaling A = ^/Tz fi = ^frw and obtain 



Kp{x,y:T) 



It I dw / dz- 



^•»^{w)-T^{)y{z) 



z — w 



i)x{z) 



The saddle points for JRi? are 



0{t 



-9/2N 



T - 2^ + ^'- 



1 fx 



(B.l) 

(B.2) 



2 Vr 



+ 0(r-9/2) 



d'Azc) = ^ ze = <! 1 - |r-i + 0(r-6/2) ^^^^^^ ^ ^ 



^r-i+0(r-«/2) 



©(r- 



-A^--,+0{r 



-9/2- 



-9/2- 



(B.3) 



The real part of 'dx{w) has local maxima at Wc — ±1 in the vertical direction, which is suitable for 
the contours ^l.r] viceversa the real part of —dy{z) has local max. at z ~ again in the vertical 
direction, suitable for iM.. The usual steepest descent asymptotics then gives. 



\ T 



(B.4) 



with C a constant (independent of x, y) whose value is immaterial for our considerations. As a 
consequence of this simple estimate we immediately deduce that 



Indet Id-X^(»,»;T) 






<^(6-a)"C"T^te-5^ = - 
— 1 r 



(h — a)CT 2c i'^ 



„^i .-{h-a)CT 2e- 

and hence the Fredholm determinant tends (exponentially) to one. 



(B.5) 
(B.6) 



C On the small— norm theorem 

We briefly review the outline of the proof of a "small norm theorem" in a somewhat idealized 
situation to sec how the length of the contours affects the estimates. We consider a simplified 
situation where the jumps arc only on the circles (which are anyway the main source of correction 
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terms). Consider the RHP for £{z) on the a coUection of circles S. For simphcity let's assume 
Y. = Sr with Sr '■= {\z\ = r}. 

£+(z)=5_(z)(l + G(z)), \z\^r (C.l) 

£{z) = l + 0{z-'^) (C.2) 

G{z) is assumed to be in any L^ of S and in fact we assume it also smooth (this is a very relaxed 
assumption but it fits our setup and most setups). The solution can be presented as 

fS^^um^ (C.3) 

J (W — ZjZlTT 

Taking the boundary value from the right we have the following singular integral equation for 

p{z) ■.^£-{z)-leL^iSr) 

{Id-C-{*-G))p = C-{G) (C.4) 

where C_() is the Cauchy boundary value operator on S. On the circle this is the projectors on the 
analytic part inside the disk, and • • G stands here for the (right) multiplication operator by the 
matrix G(z). It is easy that the operator norm in L'^{Sr) of C_ for the circle is one (and independent 
of the radius!). Therefore the operator norm of C := C_(» • G) is bounded by 

infill <||G||oo (C.5) 

Thus C is invertible as long as ||G||oo is less than one and 

iiia^-^r^iii< , ii^ii (C.6) 

Thus the solution to the singular integral equation has a bound 

||p||2 < \\\{Id-L)-'\\\ ||C_(G)||2 < ^ "^^11 (C.7) 

Now (here \M\ is any norm on the space of matrices for example the Hilbert-Schmidt norm |Mp = 
Tr(MM^) while \dw\ is the arc-length) 

_ IIGill, + ll^lbllClb ^ _^J^^ fiidl, + _HL) < (C,9) 



27rdist(z, 5^) - 27rdist(z,5r) V 1 - l|G| 

What is important to us here (this is only a rough estimate) is that the estimate is proportional to 
the length of the circle £{Sr)- In our problem these are proportional to A^"' and hence they improve 
the estimate of the error for z in the neighborhood of infinity (which is what we need). 
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